Macho Proper Motions 
From Optical/Infrared Photometry 



Andrew Gould 

Dept of Astronomy, Ohio State University, Columbus, OH 43210 

Douglas L. Welch 

Dept of Physics & Astronomy, McMaster University 
Hamilton, ON L8S 4M1 Canada 

gould@payne.mps.ohio-state.edu; welch@physics. monaster. ca 

Abstract 

Optical/infrared photometry can double the number of proper motion measure- 
ments of Massive Compact Objects (MACHOs) relative to single band photometry. 
The proper motion of a MACHO can be measured by finding the ratio q of the 
(known) radius of the source star to the Einstein radius of the MACHO, q = 6 S / ' O e . 
A classic method for doing this is to look for the effect on the light curve of the 
hnite size of the source. A modification of this method proposed by Witt (1995) is 
to look for color changes in the light curve due to the fact that the limb darkening 
of the source is different in different bands. We demonstrate that the "classical" 
method is not feasible unless the MACHO actually transits the source: if the MA- 
CHO passes at say 1.5 source radii, there is still a sizable ~ 5% effect, but the 
light curve cannot be distinguished from point-source light curves with different 
parameters. However, color measurements in V (0.55 /mi) and H (1.65 /mi) reduce 
the errors by a factor ~ 120 and permit proper motion measurements at impact 
parameters of up to 2 source radii. Color maps in V — H are also useful in the 
detection of planetary systems. Giant stars have a "red ring" in such maps. A 
planet which transits this ring gives rise to a distinctive signature which can help 
in the measurement of the planetary system's proper motion. 
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1. Introduction 



The MACHO collaboration (Alcock et al. 1995; Bennett et al. 1995) and the 
OGLE collaboration (Udalski et al. 1994) have found a total of 56 candidate Mas- 
sive Compact Object (MACHO) events toward sources in the Galactic bulge. The 
events are detected by measuring the flux of the source st function of time. 

The flux is magnified by a factor A[x(t)] (Paczyhski 1986) 

. / x I 2 , , 

AW = ^Ti)^' (L1) 

where x is the separation between the lens and source in units of the Einstein 
radius, 6 e , 



4GMB LS 



c 2 DqlDos 



M is the mass of the lens, and -Dol? -^LS? an d -^OS are the distances between 
the observer, lens, and source. If the lens is moving with speed v relative to the 
observer-source line of sight, then 



x(t) = ^(t-t ) 2 + /3 2 , (1.3) 

where to is the time of maximum magnification, f3 is the impact parameter in units 
of e , and is the Einstein ring crossing time, 

v 

The only information that one recovers from the event about the lens itself is the 
time scale, which is a complicated combination of the three physical parameters 
that one would like to know, Af, -Dol? an d v - (The distance to the source is 
assumed known from spectroscopic parallax.) In order to understand the nature 
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of these events, one would like to gain additional information. If one could measure 
the "proper motion", //, 



and the "parallax" (from which one determines the projected speed v = [Dqq/ Dis]v), 
one could obtain individual masses, distances, and speeds of the objects (Gould 
1992, 1995c). One parallax (of an unusually long event, ~ 4 months) has been 
measured from the ground (Bennett et al. 1995), and parallaxes of lensing events 
with bulge giant sources could be measured routinely using a satellite telescope 



A number of methods have been advanced for measuring proper motions. One 
method is based on the fact that if the lens transits (or nearly transits) the face of 
the source, then the light curve deviates from standard form (1.1) (Gould 1994a; 
Nemiroff & Wickramasinghe 1994; Witt & Mao 1994). Other methods have been 
advanced by Maoz & Gould (1994), Hog, Novikov, & Polnarev (1994), Simmons, 
Willis, & Newsam (1995), Loeb & Sasselov (1995), and Han, Narayanan, & Gould 



The most widely applicable of all of these methods is the photometry of transit 
events. For giant sources, the mean source radius is (r s ) ~ 22 (Gould 1995b). 
For events generated by bulge lenses, the Einstein ring (projected on the source 
plane) is ~ 400 r Q (M/0A M ) 1 / 2 (DLs/kpc) 1 / 2 . Hence of order 5% of these events 
are transits. In addition to being big, bulge giants are ideal sources because they 
are bright and hence can be easily and accurately photometered. In addition, as 
mentioned above, a future satellite could measure parallaxes for giants. 

Even if the lens does not transit the star, the light curve is still affected. As 
we will show below, the magnification is then 



fi = LO0, 




(Gould 1994b, 1995a). 



(1995). 
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where 

q= h (1 - 7) 

6 S is the angular radius of the source, and A ~ 1 is a factor that depends on the 
limb darkening of the source. At x = l.bq for example, the flux is augmented 
by ~ 5%. Gould (1994a) argued, however, that this effect would be essentially 
undetectable because it could masquerade as a change in the other parameters. 
Only if the lens actually transited the source would there be enough structure in 
the light curve to detect the finite size of the source. 

Witt (1995) pointed out that limb darkening would actually generate color 
terms in the light curve. The sense of the effect of limb-darkening is that for red- 
der wavelengths, the star will more nearly resemble a uniform disk. This implies 
that a V — H color map of a star has a bright red ring. Hence the flux increment 
in equation (1.6) is greater when the star is observed in redder bands. Since gravi- 
tational lensing is usually achromatic, this color term is an unambiguous signature 
that the lens is resolving the source. 

Here we confirm the claim by Gould (1994a) that the deviation of the light 
curve cannot be detected in a single band when f3 > q. But we also show that two- 
band (specifically V and H) photometry increases the sensitivity by a factor ~ 120, 
corresponding to an increased statistical power ~ 10 4 . Optical/infrared photome- 
try enables the detection of proper motions even for /5 ~ 2g, thus about doubling 
the number of proper motion measurements relative to single-band photometry. 
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2. Magnification of Limb-Darkened Stars 



Consider a ring of light of radius y6 e whose center is separated by x6 e from 
the lens, with i < 1. A point at <f> along the ring will be separated from the lens 
by r = (x 2 + y 2 — 2xy cos c/)) 1 ^ 2 . Since r<l, A(r) — > r _1 . Hence the ring will be 
magnified by 

A(x,y) = -^ jd<P(x 2 + y 2 -2xyco S <P)- 1 / 2 = A(x) 1 + + ^ (f ) +••• 

(2. 

Limb darkening of the surface brightness S of stars can be parameterized by 



Limb-darkening coefficients for a cool (4500 K) giant (log*/ = 1.5) of solar metal- 
licity in V and H are k\ = 0.798, t% = -0.007, «f = 0.206, and K f = 0.331 
(Manduca, Bell & Gustafsson 1977; Manduca 1979). Integrating the flux over the 
entire star [and keeping terms in equation (2.1) only to second order] yields a total 
magnification given by equation (1.6), with 

A = M-""-*", (2 .3) 

30 - 10«i — 5« 2 

Note that for j < 1 equation (2.3) remains valid even for x <; 1 (Gould 1994a). 
For V and H the limb-darkening corrections are = 0.86 and A H = 0.93. (For 
U, B, V, I, H, and K, A is respectively 0.77, 0.82, 0.86, 0.89, 0.91, 0.93, and 
0.95. In this paper we adopt V and H as the "blue" and "red" bands in order to 
obtain the largest leverage in A blue — A red while meeting practical observational 
constraints. We choose H rather than K to avoid the greater sky noise in K. We 
choose V rather than U or B, first because during bright time most bulge giants 
are fainter than the sky in the latter bands, and second because many bulge fields 
are heavily obscured. We note nevertheless that for some fields and some observing 
conditions, other combinations of colors may be better.) 
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3. Statistics of Light-Curve Fitting 



The flux from a lensed star of finite size can be written as a function of six 
parameters a n = (to, /3, lo } i*o, B } e) 

e 



F(t) = F A[x(t;t ,P,u)] 1 + 



x 



+ B, 



3.1 



where Fq is the flux from the source without lensing, B is the light from any 
unresolved and unlensed star such as the lens itself or a binary companion to the 
source, and e = Aq 2 /8. Note that B must be included as a free parameter even 
if the best fit is consistent with no unlensed sources because such sources cannot 
be ruled out a priori. For simplicity of exposition we set the zero point of time to 
to = and normalize the unit of time so that u = 1 and the unit of flux so that 
Fq = 1. We assume that the best fit value of i? <C 1, and that e <C 1. Suppose 
that a series measurements of F{t k ) are made at times with uncertainties <r&, 
and that equation (3.1) is fit to the measurements. The covariance matrix of the 
errors of this fit is given by c tJ where 

k 



b- 1 



da l daj 



(3.2) 



We evaluate c tJ for f3 = 0.1 assuming that measurements are made at a rate 
Nlo with (photon-limited) precision a = (JqA 1 ! 2 from t = to — (when the event 
is first noticed) until t = to + 5cu _1 . We find 

= iV" 1 / 2 ^_( 0.0013 3.0056 0.2683 0.0867 0.0963 0.2978 ), 



( 1 
-0.015109 

-0.012763 

-0.000661 

0.002785 

V -0.015135 



-0.015109 
1 

0.985820 
0.671057 
-0.772565 
0.999998 



-0.012763 
0.985820 
1 

0.784270 
-0.865561 
0.985477 



-0.000661 
0.671057 
0.784270 
1 

-0.987385 
0.669518 



0.002785 
-0.772565 
-0.865561 
-0.987385 
1 

-0.771250 



-0.015135 \ 
0.999998 
0.985477 
0.669518 

-0.771250 



1 

(3.3) 
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where c tJ = c tJ / yjcacjj is the matrix of correlation coefficients, and where the 
indices refer to the vector components (to, /3, lo } Fq } B } e). The measurement error 
of e is very large, corresponding to a fractional error in q 

q Ac;-' 0.01 V .J ' ' ' 

Clearly when f3 > q (the regime for which this calculation is valid), q cannot be 
measured unless the photometry is extraordinary. 

The reason for the poor precision is the extreme degeneracy of e with the other 
parameters, most notably with fi for which 1 — c# e = 2 X 10 -6 . That is, the finite- 
source light curve (e > 0) is very well mimicked by a point-source light curve with 
parameters adjusted by 

8a l = — e-^-. (3-5) 

If the a n are held fixed for i ^ 6, then the light curves with and without the finite- 
source term are easily distinguished. But, if the remaining parameters are adjusted 
by equation (3.5), the curves are essentially the same. See Figure 1. 

Let us now imagine that the observations are carried out simultaneously in 
V and H with, for example, an optical/infrared camera equipped with a dichroic 
beam splitter. For simplicity we assume that the errors are equal in the two bands. 
If the V and H light curves are initially fit separately, there will be 12 parameters, 
two for each of the parameters in the single band case. The covariance matrix 
will then be block diagonal, with each 6x6 block a duplicate of equation (3.3). 
However, we may now impose the following four constraints: 

toy = to,H, /V = Ph, u v = uj h , 0.93ey = 0.866^. (3.6) 

The last constraint arises because the finite-size term differs in the two bands by a 
known amount. As a practical matter, each constraint is imposed by successively 



7 




Figure 1. Theoretical light curves for a lensing event with finite source size q = 0.067 and [3 = 0.1. 
Shown are the true curve (solid), the best fit curve assuming a point source but with the remaining five 
parameters freely adjusted (short dashes), and the point-source curve with the remaining parameters held 
fixed (long dashes). Note that although the true curve can be easily distinguished from the point-source 
curve with fixed parameters, it is indistinguishable from the best-fit point-source curve. 



adjusting c tJ by 

, Cikai-Cjiai 



(3.7) 



where a l a l = is the constraint. For example, the last constraint is represented 
by a l = (0, 0, 0, 0, 0, 0.93, 0, 0, 0, 0, 0, —0.86). Carrying out this procedure, we then 
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find that the covariance matrix is given by 



N -i/2_ZH_ (0.0009 0.0264 0.0323 0.0462 0.0445 0.0025 ), 
0.01 v ; ' 
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(3.8) 



Note that the error in e has fallen by a factor ~ 120 relative to equation (3.3). 
That is, the "120" in equation (3.4) is replaced by unity. Thus for a lens that has 
an impact parameter of two stellar radii (/3 = 2q), q can be detected at the 3 a 
level with N = 150 and ctq = 1%. The mathematical reason for this enormous 
improvement is that while the light curve in each band can be well fit with a point- 
source curve by suitably adjusting the remaining parameters according to equation 
(3.5), these adjustments are different for V and H . Hence, the adjustments of the 
other parameters cannot simulate the change in both bands simultaneously. From 
a physical standpoint, the finite source term introduces color terms which cannot 
be mimicked by ordinary lensing of a point source, because lensing is achromatic. 

This last point deserves closer examination. In fact, ordinary lensing can have 
color terms if the source is a blend, B ^ 0. If for example the lensed source is bluer 
than the unlensed blended light, then the event will become bluer as it approaches 
the peak. However, strong constraints can be obtained on the colors of the lensed 
and unlensed sources by comparing the color of the event at x 3> 1 and at x ~ 0.5. 
Note that this comparison is not significantly affected by the finite size of the source 
because typically q ^ 0.1. Once the colors of the the lensed and unlensed sources 
are known, it is trivial to predict the expected color at the peak for a point-source 
event. Deviations from this prediction then must be due to a finite source. 
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Of course, the fitting procedure automatically takes all this into account, but 
the physical explanation is nonetheless useful. We showed that N = 150 mea- 
surements per Einstein crossing time are necessary to detect q at q = 2/3. This 
implies a total of 900 measurements during the time interval — < t < 5iu _1 . 
But according to the above analysis, not all of these measurements are equally 
necessary. To test this, we imagine that the measurements are made at N = 150iu 
for 0.5 crossing times around the peak, around x = 0.5, and around x = 4.7, and 
that during the remaining time the observations are conducted at only N = 
We find that this reduction from 900 to 290 measurements increases the errors by 
only ~ 40%. 

In our illustrative example, we chose (3 = 0.1. How does sensitivity to q depend 
on /3? To address this question, we write 



and evaluate rj for various values of (3. We find for (3 = 0.05, 0.10, 0.15, 0.20, 
and 0.3, that rj = 0.54, 0.98, 1.26, 1.43, and 1.84. That is, color measurements 
of proper motions depend mainly on the number of stellar radii of the impact 
parameter f3/q, and only weakly on (3 and q separately. 

We conclude that measuring proper motions of lensing events from optical/infrared 
colors is substantially more effective than similar measurements in a single band, 
which latter are possible only for transits. The rate at which proper motions can 
be conclusively detected using a combination of optical and infrared photometry 
is at least twice the single bandpass rate. 




(3.9) 
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Figure 2. Change in color as a function of source-lens separation in units of the stellar radius, assuming 
e ^> 8 S - The zero point of the color is that of the unlensed source. 



4. Transits and Planetary System Lenses 

If the lens transits the face of the star, then the light curve in each band 
deviates substantially from the limiting form of equation (1.6). In this case it is 
straight forward to measure the proper motion from each band separately. Colors 
are not needed. However, the color term does help confirm the transit. We show 
this color (relative to the color of the unlensed star) as a function of source-lens 
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separation in Figure 2. 

The main value of colors during transits is in the search for planetary systems. 
If a star with a small planet having mass m <C M acts as a lens, then for most 
of the lensing event, the light curve will look like a standard microlensing light 
curve. However, as first pointed out by Mao & Paczyhski (1991) in some cases 
the planet can dramatically alter the light curve during a fraction (m/Af) 1 / 2 of 
the Einstein crossing time. Gould & Loeb (1992) analyzed the problem in detail 
but restricted their attention to the case of point sources. They estimated for 
planets in the "lensing zone" [with projected separations from the parent star 

3-6 AU (M/M©) 1 / 2 ] that the planet would give rise to a perturbation <; 5% 
in ~ 17%(m/mj up iter) 1 ^ 2 of the events. Note that the Einstein radius associated 
with a Jupiter mass planet is ~ 50 tq / Rq where is the radius of the Sun and 
i?o ~ 8 kpc is the galactocentric distance. Hence the point-source approximation 
remains valid even for bulge giant sources with radius r ~ 10 tq. 

However, for smaller planets, particularly Earth mass planets, the point-source 
approximation breaks down. In this case, one can detect the planet only if the 
caustic associated with the planet actually transits the star. Even then the degree 
of brightening may be quite modest. One can show that if the planet is separated 
from the star by more than one Einstein radius, then the star brightens by a factor 
A ~ 1 + 2(m/M)(9 e /9 s ) 2 . For an Earth -mass planet at 2 kpc from the Galactic 
center, this is A ~ 1 + O.O2(r/2Or0) 2 . The effect is therefore detectable even for 
large sources. This is important because the probability that the planet will transit 
the star ~ (r / Rq) / Q e . The mean radius of giants (Mj < 0.5) in the Galactic bulge 
is (r) ~ 22 (Gould 1995b). The probability that a planet will transit such 
a star is ~ 2%(M/Mq) 1 / 2 . Given that an aggressive search could detect ^ 100 
giant-source events per year, there is a significant chance to detect some small 
planets. 

If such planets are detected, one would like to know as much about them as 
possible. The peak excess light A — l constrains the quantity (m/M)(9 e /9 s ) 2 . If one 
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Figure 3. Change in V — H color of the excess light (A — 1) due to a very small lens passing over the 
face of the star, 6 e <C 9 S . The zero point of color is that of the unlensed star. Note that there is a large 
color shift for 0.9 0/6 s < 1 which gives rise to a "red ring" in the color map. This ring is helpful in the 
analysis of events generated by planetary systems with small (e.g. Earth-mass) planets. 



knew the geometry of the transit (i.e., the impact parameter in units of the source 
radius) then the duration of the transit would yield (6 s /6 e ) and hence the proper 
motion // = u6 e = v/Dq^. The question is, how does one know the geometry of 
the transit? In principal, this can be extracted from single band photometry by 
measuring the length of time it takes the caustic to move across the edge of the star, 
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that is the time from no excess magnification to full excess. This [together with 
the above mentioned measurement of (m / M) 1 ! 2 9 e ] gives the angle between the 
edge of the star and the trajectory of the lens. In practice, much more information 
would be available if the small lens were transiting a narrow ring rather than the 
edge of a disk. The enhanced sensitivity of a ring has been pointed out by Loeb & 
Sasselov (1995) who showed that such a ring exists in narrow band Call images. 
As shown in Figure 3, V — H color maps also have a narrow ring. As with Call 
images, this ring can help resolve the proper motion of planetary systems. Under 
most observing conditions V — H color maps will have substantially higher signal to 
noise than Call images in part because they are based on broad-band photometry 
and in part because V band is less heavily extincted than the 393 nm Call line. 
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